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ABSTRACT 


The  model  of  Haskell  for  explosion  source  time 
functions  and  spectra  fails  to  satisfy  data  in  the 
short-period  band  recorded  teleseismically  from  the 
three  Amchitka  Island  underground  nuclear  tests: 

LONG  SHOT,  MILROW  and  CANNIKAN.  A  more  recent  model 
due  to  Mueller  and  Murphy  satisfies  the  data  quite 
well.  The  difference  in  the  two  models  is  basically  in 
the  falloff  at  high  frequencies.  A  simple  revision  of 
Haskell's  model  produces  waveforms  and  spectra  nearly 
identical  to  ones  from  Mueller  and  Murphy's  model. 

This  revision  requires  velocity  waveforms  to  have  a 
rise  time  of  extremely  short  duration  at  the  elastic 
boundary,  a  premise  validated  by  actual  near-field 
measurements. 

Waveforms  are  derived  from  the  revised  Haskell 
model  and  the  Mueller  and  Murphy  mod^l  and  illustrated 
for  pressure  at  the  elastic  boundary,  reduced  displace¬ 
ment  potential  at  the  elastic  boundary,  and  far-field 
displacement.  Corresponding  spectra  are  derived  and 
illustrated. 
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SYNOPSIS  OF  PREVIOUS  WORK  ON  SOURCE 
SPECTRA  AND  YIELD  SCALING 


The  response  of  an  infinite  homogeneous  elastic 
medium  to  a  pressure  function  acting  on  the  boundary  of 
a  spherical  cavity  has  a  well-known  solution  (e.g,, 
Sharpe,  1942;  Blake,  1952).  The  solution  for  simple 
time  functions  of  pressure  such  as  a  step  or  a  decay¬ 
ing  pulse  are  straightforward.  Recently  attempts  have 
been  made  to  model  the  exact  source  time  function  for 
nuclear  detonations  with  the  aid  of  close-in  empirical 
measurements.  Among  other  reasons,  we  feel  it  is 
imperative  to  obtain  a  nearly  exact  model  for  the  source 
function  so  that  spectral  ratios  applied  to  the  short- 
period  band,  from  0.2  to  5  cps,  can  be  formed  and 
analyzed  in  an  intelligent  manner  for  purposes  of 
discrimination. 

Haskell  (1967)  formulated  the  source  spectrum  of 
an  explosion  by  fitting  a  parameterized  function  to 
actual  displacement  potential  functions  calculated  from 
measurements  just  outside  the  elastic  radius  of  under¬ 
ground  nuclear  explosions.  The  function  used  results 
in  a  far-field  displacement  spectrum  which  is  asymptotic 
as  the  inverse  of  the  fourth  power  of  frequency  (von 
Seggern  and  Lambert,  1970).  This  falloff  entails  a 
high-frequency  scaling  of  displacement  inversely  pro¬ 
portional  to  the  cube  root  of  yield  (i.e.,  large 
explosions  emit  less  high  frequency  energy  than  small 
explosions)  and  a  low-frequency  scaling  of  displacement 
proportional  to  yield. 
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Mueller  and  Murphy  (1971)  have  formulated  a  nodel 
baaed  on  theoretical  considerations  of  the  nediun 
response  to  an  explosive  source  near  the  surface.  The 
theoretical  underpinnings  are  basically  those  of  Sharpe 
(1942)  and  Blake  (19S2).  The  far-field  displacement 
spectrum  from  these  models  is  asymptotic  as  the  inverse 
of  the  square  of  frequency  for  high  frequencies.  This 
entails  a  high-frequency  scaling  of  displacement 
directly  proportional  to  the  cube  root  of  yield,  in 
contrast  to  llaskell*s  model.  The  low-frequency  scaling 
is  the  same  as  for  Haskell's; that  is,  displacement  is 
proportional  to  yield. 

The  purpose  of  this  report  is  to  compare  the  two 
models  to  data  taken  from  teleseismic  recordings  of 
the  three  Amchitka  Island  underground  tests.  A  full 
development  of  explosion  functions  will  be  made  after 
reviewing  the  data;  this  will  include  representations 
for  the  pressure  and  displacement  potential  at  the 
clastic  radius  of  the  explosion,  the  far-field  dis¬ 
placement,  and  the  spectrum  of  each  of  these  functions. 
In  this  report  we  exclude  the  effects  of  the  free 
surface  on  the  theoretical  source  functions,  and  for 
comparison  with  theory  we  attempt  to  remove  the  free 
surface  effects  in  the  data. 


BODY-NAVE  SCALING  BOR  AMCHITKA  ISLAND  TESTS 


Scaling  of  flrat  notion  aaolitudes  for  LONG  SHOT. 

MILRON.  and  CANNIKAN 

Na  will  conpara  tha  aaasurad  ralativa  aaplitudes 
of  tha  P-wava  for  tha  thraa  Anchitka  Island  tests  with 
values  calculated  froa  tha  two  theoretical  scaling 
foraulas  under  consideration.  Tha  reported  yields  of 
the  three  explosions  are:  LONG  SHOT,  SO  kt;  MILRON, 

1000  kt;  CANNIKAN,  S000  kt.  In  order  to  avoid  the 
effects  of  the  reflection  fron  the  surface,  which  is 
delayed  approxinately  a  half  second  relative  to  the 
initial  pulse  for  LONG  SHOT  (Cohen,  1969)  and  nore  for 
the  other  two  shots,  we  aust  use  the  unbiased  anplitudes 
of  the  initial  upward  ground  notion  at  the  receiver. 

By  neasuring  this  first  quarter- cycle  at  several  common 
stations,  von  Seggern  and  Laabert  (1971)  deterained  the 
ratio  of  MILRON  to  LONG  SHOT  aaplitude  to  be  6.49.  Ne 
have  repeated  the  sane  procedure  with  CANNIKAN  and 
MILRON  at  seven  covion  stations  for  which  data  could 
be  obtained  at  the  Seisaic  Data  Laboratory;  the  aeasure- 
aents  are  listed  in  Table  I.  The  average  ratio  of 
CANNIKAN  to  MILRON  aaplitude  is  thus  estiaated  to  be 
2.56. 

Froa  Haskell's  scaling  theory  for  granite,  we  can 
calculate  using  the  displaceaent  spectrua  formula  as 
given  by  von  Seggern  and  Laabert  (1970)  the  relative 
anplitudes  at  a  frequency  of  one  cps.  Froa  Haskell's 
work,  one  is  not  able  to  account  for  the  effects  of 
snail  changes  in  aediua  properties.  The  granite  scaling 
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TABLE  I 


Cannikan 


Station 

Miirow 

LAO 

2.59 

KN 

2.36 

RK 

2.63 

CR 

2.47 

HN 

2.01 

BE 

3.16 

TFO 

2.70 

Average 

2.56 

Standard  Deviation 

.35 

35 


is  used  because  the  Amchitka  tests  apparently  follow  the 
"hard  rock"  magnitude-yield  scaling  at  NTS  (von  Seggern 
and  Lambert,  1971),  The  amplitude  ratios  arc  calculated 
at  one  cps  because  this  agrees  fairly  well  with  the 
period  of  the  first  cycle  of  motion  at  all  the  stations 
for  all  three  shots.  There  is  a  slight  shift  to  lower 
frequencies  for  increasing  yield  though.  The  scaling 
ratio  between  two  shots  in  granite  at  1  cps  is  given 
by: 


u?  Y2  |1  ♦  4tt2G2Y§/3  A211/2  (1  ♦  4tt2(;2yJ/3j  5/2 
U1  *T  (1  HuVy^  A2]1^  (1  ♦  4tt^G2y|^3) S/2 


where 

u  ■  vertical  displacement  in  far  field 
Y  ■  yield 

G  *  .0185  -  constant  for  granite 
A  ■  1  ♦  24B 

B  *  ,24  -  constant  for  granite. 

Calculations  using  (lj  give  a  ratio  of  MILROW  to  LONG 
SHOT  amplitude  of  6.58  and  a  ratio  of  CANNIKAN  to 
MILROW  amplitude  of  1.51, 

The  scaling  theory  of  Mueller  and  Murphy  (1971) 
involves  considerably  more  terms  which  must  be  esti¬ 
mated  from  available  data.  Their  scaling  relation 
between  two  shots  in  different  media  at  different 
depths  at  1  cps  is  given  by: 


1/2 


u. 

j 

U, 


rcl2  C2^ 
reli  clu2 


a?  ♦  4tt^  4it 
’  T”1^  *  - 

<*2  ♦  j  4-tt 


1  r«5  * 


I 


+  1 


picj 

fZ 


f 
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i/: 


where 
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c 

M 

W< 

e 

ct 

P 


os 


Koc 


j  16it^6|  ♦  (l*2Bj)4t^Qj  ♦  wqj 

I  16*^2  ♦  (1-282)4ttZwo2  4  w02 


clastic  radius 

compress ional -wave  velocity 

rigidity  modulus 

c/rei 

(X  ♦  2u) /4m  (X  is  Lame  *  s  constant) 

ku>  (k  being  approximately  2  for  granite) 

1.5  pgh  (p  is  density,  g  is  gravitational 

acceleration,  and  h  is  depth) 
(4M/3)(rc/rel)3(rc  is  final  cavity  radius). 

is  the 


oc 


The  value  pQS  is  the  initial  pressure,  and  p  ^ 
residual  pressure  at  large  t.  Several  assumptions  must 
be  made  now.  We  first  assume  X  s  2p  thus  making  3=1. 
This  implies  that  Poisson’s  ratio  is  equal  to  .33,  a 
value  more  suited  to  rocks  at  shallow  depth  than  .25. 
We  assume  that  rfil  scales  exactly  as  the  cube  root  of 
yield,  and  that  rc/rel  equals  0.1.  Both  these  assump¬ 
tions  are  reasonable  approximations  (Mueller  and 
Murphy,  1971).  For  the  LONG  SHOT  site,  we  use  c  =  3.5 
km/sec  as  indicated  by  the  velocity  log  at  the  LONG 
SHOT  site  (Lambert  et  al.,  1969).  For  the  MILROW  site, 
we  use  c  =  4.0  km/sec  as  indicated  by  the  velocity  log 
prepared  by  Snyder  (1969).  We  do  not  as  yet  know  the 


(2) 


H 


the  velocity  at  the  CANNIKAN  site;  but,  using  an  accepted 
velocity-depth  law  of  v  ■  vQ  ♦  a(h-hQ)  where  v  is  the 
estimated  Velocity  for  the  CANNIKAN  test  at  1,8  km  depth, 
vQ  is  the  MILROW  velocity  at  its  1.2  km  depth,  and  - 
a  ■  1  km/sec/km,  we  obtain  a  velocity  of  4.6  km/sec  for 
CANNIKAN.  Using  the  relation  X  ■  2y,  we  can  calculate 
the  rigidity  moduli  for  all  three  sites:  LONG  SHOT,  7.4; 
MILROW,  1.6;  CANNIKAN,  12,7.  To  calculate  the  peak  pres¬ 
sures,  pQS,  we  use  p  -  2.4  gm/cc  at  all  three  sites  as 
indicated  by  the  fairly  constant  density-depth  profile 
at  the  LONG  SHOT  site  (Lambert  et  al.,  1969).  The 

remainder  of  the  parameters  in  (2),  a,  w  ,  and  p  ,  can 

0  oc 

now  be  calculated.  Application  of  (2)  then  to  the  three 
tests  results  in  predicted  amplitude  ratios  of  MILROW 
to  LONG  SHOT  of  7.57  and  of  CANNIKAN  to  MILROW  of  2.36. 

Table  II  summarizes  the  observed  and  calculated 
amplitude  ratios.  Both  the  scaling  relations  of  Haskell 
and  Mueller  and  Murphy  agree  with  the  actual  MILROW  to 
LONG  SHOT  amplitudes  quite  well.  For  CANNIKAN  to  MILROW, 
however,  Haskell's  scaling  is  significantly  in  error 
while  that  of  Mueller  and  Murphy  agrees  very  well  with 
the  data.  Together  with  the  RK-ON  spectral  ratios  for 
LONG  SHOT/MILROW  due  to  Sax  (1972),  it  was  in  fact  the 
failure  of  Haskell's  scaling  theory  to  correctly  predict 
the  CANNIKAN  amplitude  which  initiated  this  study. 

For  three  reasons  the  above  comparisons  cannot  be 
entirely  convincing  in  supporting  Mueller  and  Murphy's 
scaling  relation  over  Haskell's.  Firstly,  the  change  in 
the  spectral  shape  with  increasing  yield  in  conjunction 
with  the  frequency  response  of  the  recording  systems 


TABLE  II 


Mil  row  Cannikan 

Long  fihot  Mi  lrow 


Ratio 

Log  (Ratio) 

Ratio 

Log  (Ratio) 

Measured  Data 

6.49 

.81 

2.56 

.41 

Haskell's  Model 

5.61 

.75 

1.28 

.11 

Mueller  5  Murphy1 
Model 

f  s 

7.57 

.88 

2 . 36 

.37 

may  distort  the  measured  relative  amplitudes  at  1  cps; 
however  we  have  estimated  that  this  effect  cannot  be 
more  than  0.1  or  0.2  in  the  logarithm.  Secondly,  the 
number  of  assumptions  and  estimations  required  in 
applying  Mueller  and  Murphy's  scaling  relation  calls 
for  considerable  error  limits  on  its  results.  Thirdly, 
Haskell's  theory  allows  for  no  comparison  between  shots 
at  different  depths  and  in  media  other  than  the  four 
which  he  examined  and  is  therefore  not  a  general  scaling 
theory. 

P-wavc  spectral  ratios  at  RK-ON 

To  reinforce  the  results  of  the  comparison  using 
initial  P-wavc  amplitude  data,  we  will  examine  the 
entire  short  period  hand  of  the  recorded  P-wavcs  from 
MILROb  and  LONG  ‘•‘HOT  and  compare  this  with  spectral 
calculations  from  the  two  scaling  relations. 

Specifically,  we  examine  the  spectra  of  LONG  SIKIf 
and  NILROh  P  waves  at  ML -OS  as  shown  hy  von  Seggern  and 
Lambert  (1971).  These  unsmoothed  spectra  are  corrected  for 
noise  hy  subtracting  out  the  spectrum  of  a  noise  sample 
just  preceding  the  signals.  RR-0N  was  chosen  among  the 
many  stations  in  that  report  because  the  signal  Is  above 
the  n«ise  out  to  S  cps  at  least  for  LONG  SIW1,  and 
because  the  modulation  of  the  spectra  by  the  surface 
reflection  is  nearly  ideally  shown  in  both  cases,  be 
can  remove  the  effect  of  the  surface  reflection  when  we 
merely  smooth  over  these  modulations  hy  connecting 
adjacent  spectral  peals  with  straight  lin^s  to  get  the 
spectra  shown  In  figure  I.  from  these  smoothed  spectra. 


we  calculate  amplitude  ratio*  of  MIMtOM  to  LONG  SHOT  at 
O.s  cp«  Interval*  and  plot  then  a*  In  Figure  2.  In 
addition,  the  amplitude  ratio  at  long  periods  can  ho 
obtained  froia  the  natch  filter  aaethod  for  Mayloigh  waves s 
von  t/eggern  and  Lambert  report  this  to  be  S.S  at  teK*ON, 
and  we  plot  this  at  .01  cps  In  Figure  2.  Using  llaskell** 
sealing  relation  for  Id  and  IttOi)  kt  shots  In  granite, 
we  compute  the  spectral  ratios  shown  in  Figure  2*  Using 
a  revision  of  Haskell**  model  to  give  a  falloff  at  high 
frequencies  proportional  to  the  Inverse  of  frequency 
squared  as  In  Mueller  and  Murphy's  model,  we  calculate 
the  ratios  shown  in  Figure  2  which  agree  remarkably 
well  with  the  actual  data*  (Ibis  revision  will  be 
discussed  in  the  nest  section,  but  Its  form  is  nearly 
equivalent  to  Mueller  and  Murphy**,  and  a  full  eaten* 
sion  of  equation  (2)  across  the  band  from  0.$  to  S  cps 
would  produce  a  line  nearly  identical  to  the  solid 
one  shown  In  Figure  2.) 

the  results  of  comparing  spectral  ratios  at  one 
station  and  first  motion  amplitudes  at  several  stations 
together  affirm  the  validity  of  the  model  of  Moeller 
and  Murphy  and  cast  doubt  on  the  validity  of  Haskell's 
model.  In  the  neat  section  we  will  delineate  toe 
essential  physical  differences  of  the  eaplosive  source 
mechanism  as  inferred  by  the  two  models  and  develop 
general  forms  for  representing  the  source  time  fund  ions 
and  spectra. 
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TINR  RHCTIOhS  m  SFfeCTIU 

liaiMUi  EEiiiati  agRtiiim 

llaskell  (19*7)  formulated  kis  Model  by  fitting  curves 
of  tke  typo 


tSi 


there  #H  Is  tkt  asymptotic  value  for  largo  t.  to  reduced 
displacement  potentials  calculated  from  data  taken  Just 
outside  Ike  elastic  radius  of  several  underground  nuclear 
detonations.  As  stated  ky  kin,  tke  fora  of  tkis  function 
vis  ckosen  so  as  to  Make  displacement,  velocity ,  and 
acceleration  functions  continuous  at  t  *  •  It  is  tke 
retarded  tine  referred  to  tke  elastic  boundary).  Ibis 
entails,  as  skown  by  von  Seggern  and  kanhert  t 197ft), 
tke  dependence  at  kigk  frequencies  for  tke  far- 
field  displacement  spectrin.  Haskell's  curves  according 
to  (3)  are  apparently  excellent  fits  to  tbe  data  he 
shoes,  tfe  reproduce  bis  example  for  a  granite  medium  In 
Figure  1.  However,  note  that  a;  the  beginning  of  tbe 
wavefora,  which  is  tbe  critical  aica  in  relation  to  the 
falloff  of  the  high-frequency  portion  of  the  spectrum, 
bis  analytic  curve  has  a  slope  significantly  less  than 
tbe  real  slope,  be  feel  that  Haskell's  requirement  of 
continuous  acceleration  and  velocity  at  the  elastic 
radius  is  physically  too  strong,  and  therefore  we  can 
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safoly  assume  tk«  nore  idealised  cast  where  acceleration 
•nd  even  velocity  are  discontinuous  at  the  elastic 
boundary  of  the  nodiun  surrounding  the  explosion.  These 
■ay  not  be  the  exact  discontinuities ,  but  the  junp  should 
occur  in  a  tine  span  which  is  nearly  instantaneous  relative 
to  the  frequency  band  under  consideration.  We  show  in 
Figure  4  a  reproduction  of  a  particle  velocity  aeasure- 
aont  outside  the  elastic  boundary  of  the  GASBUGGY  shot 
(Ferret,  1W0).  The  velocity  takes  a  largo  initial 
Junp  in  approximately  .01  seconds;  this  inplies  that 
for  frequencies  less  than  approximately  100  cps,  the 
spectrum  is  that  of  an  impulse,  i.e.,  flat.  In  contrast, 
Haskell  would  require  particle  velocity  to  bo  a  raap- 
llke  function  which  inplies  an  falloff  for  high 
frequencies  in  the  spectrum. 

initial  of  Tfri 

By  removing  the  quartlc  and  cubic  terms  from  (3), 
we  remove  the  constraints  of  continuous  acceleration 
and  velocity  in  the  near  field.  Then  we  propose  a 
source  model  by  fitting  the  function 


{(if  •  I  -e*kt||  ♦  kt  •  fe(kt)2) 


(«) 


to  the  same  data  as  used  by  Haskell.  For  the  granite 
case  In  particular,  we  get  values  of  k  •  1*.0  and 
•  •  2.04  by  requiring  that  (4)  fit  the  paak  of  the 
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measui  od  potential  and  the  asymptotic  long-term  value  ^  («>) . 
Several  calculated  points  of  this  curve  are  shown  in  figure  3. 

The  reduced  displacement  potential  (4)  and  the  pres¬ 
sure  function  on  the  elastic  radius  are  related  by 
(Rodean,  equation  4.17): 
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'cl 


(5) 


Using  (4)  in  (5),  we  obtain 


4pb  2'ip  (to) 


a(t) 


3  e“kt[B(kt)2  -  Akt  ♦  2B]  ♦  1 


(6) 


where  A  =  2B  +  1,  and  we  have  assumed  k  =  c/rgj 
and,  with  X  -  2 p,  k  =  2b/rel«  By  equation  4.22  of 
Rodean  (1971)  ,  we  have  for  w  small  in  the  frequency 
domain  or  t  large  in  the  time  domain: 


*(»)  - 

4pb^ 

Applying  this  relation  to  (6)  results  in: 

°(t)  -  c'kt[B(kt)2  -  Akt  ♦  2B ]  ♦  1 

-M- 


(7) 
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if  We  c°ntinue  tBe  previous  assumption  of  X  ■  2u  and  thus 
c  =  4b  .  We  desire  an  expression  for  the  far-field  dis¬ 
placement  also;  this  can  be  derived  from  the  reduced 
displacement  potential  (4)  if  we  employ  this  relation 
(Rodean,  equation  4.16): 


For  the  far-field  displacement,  we  ignore  the  second 
term  in  (9)  and  thus  have 


(10) 


Substituting  (4)  in  (10)  then  gives 


yT”T =  kt  e  "  Bktl 


(11) 


where  again  k  =  c/rej  is  assumed. 

Equations  (8),  (4),  and  (11)  give  the  pressure, 
reduced  displacement  potential,  and  far-field  displace¬ 
ment  waveforms.  The  parameter  k  is  medium  dependent;  and 
we  assume  as  did  Haskell,  that  it  scales  as  the  inverse  cube 
root  of  the  yield  Y,  B  is  a  medium-dependent;  parameter 
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which  defines  the  aaount  of  overshoot  in  reduced  dis* 
placenent  potential  wsvefom.  In  Figure  5a,  5b,  and  5c, 
we  show  (8),  (4),  and  (11),  respectively  as  a  function 
of  the  diaensionless  paraaeter  kt  for  five  values  of 
B:  0,  1,  2,  3,  and  4.  The  case  B  •  2  for  the  reduced 
potentials  in  Figure  Sb  closely  corresponds  with  the 
fit  of  (4)  to  the  granite  data  of  Haskell  in  Figure  3 
when  B  was  calculated  to  be  2.04*  it  appears  that 
equation  (4)  then  can  fit  observed  potentials  as  well 
as  Haskell's  original  function,  equation  (3).  Mot  only 
does  the  eliaination  of  quartic  and  cubic  terns  froa 
Haskell's  function  produce  a  aore  physically  satisfying 
function  as  explained  above,  but  it  does  so  witkout 
suffering  any  loss  of  fitting  capability. 

Spectra  of  the  above  three  functions  for  pressure 
at  the  elastic  radius,  reduced  potential  at  this  radius, 
and  far*  field  displaceaent  can  be  derived.  We  start 
with  the  reduced  potential  (4)}  the  Fourier  transfora 
is  straightforward,  but  tedious,  and  the  final  fora 
only  is  given: 


rrfisr  I*MI  ■  (-/k)*1  UlMMj  ♦ 
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(12) 


Rode an  (equation  4,19)  gives  the  relation  between  the 
pressure  spectrua  o(w)  end  *(w)  as 
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Figure  St*  Pressure  wave fora*  at  the  elastic  boundary  calculated  fro* 
revised  Haskell's  Model* 


Figure  Sb.  Reduced  displacement  potential  waveforms  at  the  elastic 
boundary  calculated  fron  revised  Haskell's  model. 
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The  M*a«t  of  tills  npniiiM  it 
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Substituting  (12)  U  (14)  Mi 
the  rolotlon  (7)#  wo  taro 
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tftls  using  k  •  c/rtl 

(*afa/k)2  «  nV1 
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Tko  for- fl#W  dlsplocMtnt  spoctro  cm  bo  dorlvod  with 
tho  old  of 


(IS) 


(16) 


which  Is  tho  Modulus  of  tho  Pourlor  tronsforn  of  (10). 
Using  (16)  In  (12),  wo  obtoln 


If  wo  ogoin  uso  k  ■  c/r#l.  Plots  of  (12),  (15),  ond  (17) 
oro  g Ivon  In  Figures  6o,  6b  ond  6c,  respectively,  os 
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Figure  6b.  Reduced  displacement  potential  spectra 
boundary  calculated  froa  revised  l!askeil*s  nodel. 


at  the  clastic 
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Figure  6c,  l;ar-field  displacement  spectra  calculated  from 
revised  Haskell's  model. 
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functions  of  the  dimensionless  parametor  w/k  for  fi  values 
of  0,  1,  2,  3,  and  4,  Note  that  the  required  falloff 
of  the  displacement  spectra  at  high  frequencies  is  present 
and  that  the  overshoot  in  the  displacement  spectra  is 
proportional  to  B. 

Mueller  and  Murphy’s  formulation 

Mueller  and  Murphy  (1971)  take  a  slightly  different 
approach  in  deriving  scaling  relations.  They  assume  the 
arbitrary  form 


o(t) 


oc 


_  -aw  t  ^ 
p  e  o+l 


(18) 


for  the  pressure  at  the  elastic  radius  where  p  =  P0/P0C 
and  p„  =  p„e  -  p  (p„0  and  p„„  are  defined  as  for 
equation  (2))  whereas  Haskell  began  with  an  arbitrary 
form  for  reduced  displacement  potential.  To  derive  the 
reduced  displacement  potential  from  the  pressure  func¬ 
tion,  we  transform  (18)  into  the  frequency  domain, 
apply  equation  (13),  and  inverse  transform  back  to  the 
time  domain  to  get  the  reduced  potential  function. 
Inverse  transforms  #,011  and  #,101  in  Nixon  (1960)  are 
employed  with  the  damping  coefficient  set  to  one-half. 
The  reduction  is  straightforward,  but  again  tedious, 
and  we  give  only  the  final  result: 

Mt)  =  c-wot/2  {.  577  [Cl2-a-~4)  -  1]  s  in  (,  866w  t) 

1  1-a+a 

-  [ — J  +  1]  cos  (.  866w  t)  } 
l-a+a4  0 

A  P  _-aw„t 
+  — e  0  +  1 
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The  far- field  displacement  is  obtained  by  use  of  (10)  as 
be fore: 


The  pressure  (18)  could  be  plotted  as  a  function  of  the 
dimensionless  par aiiot erau>0t,  hut  the  reduced  displacement 
potential  (19)  and  the  far- field  displacement  (20)  cannot. 
Thus  in  fitting  observed  reduced  potentials,  Mueller  and 
Murphy's  formulation  requires  the  additional  parameter  a 
to  be  estimated.  The  parameter  corresponds  to  k  in 
Haskell's  formulation,  and  the  parameter  p  here  corresponds 
to  his  2B.  In  Figures  7a,  7b  and  7c,  we  show  the  plots 
of  (18),  (19),  and  (20),  respectively,  as  a  function  of 
w0t  for  a  «  2  when  p  is  assigned  values  of  0,  1,  2,  3, 
and  4.  Several  other  values  of  a  were  used,  but  a  »  2 
is  the  value  Mueller  and  Murphy  suggest  for  rhyolite, 
which  is  compositionally  close  to  granite.  For  a  ■  2  and 
p  =  4  in  the  reduced  potentials  of  Figure  7b,  the  curve 
closely  corresponds  to  the  revised  Haskell  potential  in 
Figure  5b  when  B  is  taken  to  be  2,  approximately  the 
value  we  obtained  when  fitting  (4)  to  the  granite  data 
shown  in  Figure  3,  This  shows  the  analogy  of  p  to  2B, 

Spectra  of  the  above  functions  for  Mueller  and 
Murphy's  model  can  be  obtained  by  starting  with  the 
easily  transformable  pressure  function  (18).  The 
resulting  pressure  spectrum  is: 
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Figure  7a.  Pressure  waveforms  at  the  clastic  boundary 
calculated  from  Mueller  and  Murphy's  model. 
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Figure  7c.  Far-field  displacement  waveforms  at  the  elastic 
boundary  calculated  from  Mueller  and  Murphy’s  model. 
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The  reduced  potential  spectrum  follows  by  applica¬ 
tion  of  (14) :  1/2 


|Hw)  I  =  (<*>/%,)  ^ 
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The  far-field  displacement  spectrum  follows  from 
(22)  by  application  of  (16): 
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Note  that,  as  for  the  corresponding  time  functions, 
the  pressure  spectrum  (21)  could  be  manipulated  to  be 
a  function  of  the  variable  aw/oo  whereas  the  reduced 
potential  spectrum  (22)  and  displacement  spectrum  (23) 
cannot.  In  Figures  8a,  8b,  and  8c  we  plot  (21),  (22), 
and  (23),  respectively,  for  a  =  2  with  p  set  from  0  to 
4  again.  Comparison  of  Figure  8c  with  p  =  4  to  Figure  6c 
with  B  =  2  shows  the  similarity  of  the  revised  Haskell 
model  and  the  Mueller  and  Murphy  model  for  a  granite- 
rhyolite  medium. 
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Figure  8a.  Pressure  spectra  at  the  elastic  boundary  calculated 
from  Mueller  and  Murphy's  model. 


DISCUSSION 


Using  teleseismic  data,  we  have  shown  that  a  Model 
for  the  explosion  source  function  which  entails  a  far- 
field  displacement  spectrum  that  is  inversely  propor* 
tional  to  u  at  high  frequencies  is  the  proper  choice* 
Mueller  and  Murphy  (1971)  have  already  verified  this 
model  using  near-field  data.  Haskell's  requirement 
that  acceleration  and  velocity  to  be  continuous  at  the 
elastic  boundary  around  the  explosion  is  unnecessary; 
they  certainly  can  be  considered  as  discontinuous  for 
any  part  of  the  spectrum  capable  of  being  measured 
teleseismically.  We  mention  that  the  u  2  dependence  at 
high  frequencies  is  characteristic  of  earthquake 
models  also  (Aki,  1967 ;  Brune,  1970). 

_2  We  have  revised  Haskell's  formulation  to  obtain  an 
(o  model  and  found  that  the  fit  to  observed  potentials 
is  apparently  as  good  as  with  his  original  a)"-*  model. 

We  have  used  Mueller  and  Murphy's  formulation  for  the 
pressure  function  at  the  elastic  radius  to  derive 
reduced  potential  and  far-field  displacement  waveforms. 
Waveforms  and  spectra  for  the  two  uT2  models  are  similar. 
Observed  waveforms  can  be  fit  by  adjusting  the  para¬ 
meters  in  either  model,  and  this  will  be  the  subject 
of  a  future  report.  Especially  important  is  estimation 
of  the  parameter  p  (or  B)  which  is  the  ratio  of  the 
overshoot  pressure  pQ  to  the  residual  pressure  p  at 
the  elastic  boundary.  This  value  controls  the  peak  in 
the  far-field  spectra  and  affects  spectral  ratios  in 
the  band  around  1  cps.  These  pressures  can  be  estimated 


using  sha  relations  in  aquation  (2)  as  takan  fro*  Muallar 
and  Murphy*  Tha  paranatar  also  controls  tha  spactral 
shapa  noar  1  cps,  but  its  affact  is  pradictabla  sinca 
mo  ■  c/rfl  and  rfl  scalas  as  tha  cuba  root  of  yiald 
whan  datonation  dapths  ara  aqual  (Muallar  and  Murphy , 
1971). 
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